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ABSTRACT 

We report on a deep radio search for a pulsar wind nebula associated with the radio-quiet neutron star RX J0822- 
4300 in the supernova remnant Puppis A. The well-determined properties of Puppis A allow us to constrain the 
size of any nebula to less than 30 arcsec; however we find no evidence for such a source on any spatial scale up 
to 30 arcmin. These non-detections result in an upper limit on the radio luminosity of any pulsar-powered nebula 
which is three orders of magnitude below what would be expected if RX J0822^300 was an energetic young 
radio pulsar beaming away from us, and cast doubt on a recent claim of X-ray pulsations from this source. The 
lack of a radio nebula leads us to conclude that RX 10822^-300 has properties very different from most young 
radio pulsars, and that it represents a distinct population which may be as numerous, or even more so, than radio 
pulsars. 

Subject headings: ISM: individual (Puppis A) - ISM: supernova remnants - pulsars: general - radio continuum: 
ISM - stars: neutron: individual (RX J0822-4300) 



1. INTRODUCTION 

While the vast majority of neutron stars so far discovered 
are seen as radio pulsars, there are also a small but increas- 
ing number of neutron stars which have very different obser- 
vational properties. Approximately half these sources are soft 
7-ray repeaters (SGRs) or anomalous X-ray pulsars (AXPs), 
both of which show pu lsed X-rays at long periods (P ~ 10 s) 
(e.g. Mereghetti 1999 ). The remaining sources are grouped 
together as "radio-quiet neutron stars" (RQNS; Cara veo, Big- 
nami, & Trumper 1996; |Brazier & Johnston 1999 ), most of 
which are characterized by unpulsed thermal X-ray emission 
at a temperature of a few million degrees, a complete lack of 
radio emission, and very high X-ray to optical ratios. Many 
of these sources have been associated with supernova remnants 
(SNRs), and are thus probably quite young (<;20 kyr) objects. 

The AXPs and SGRs are quite distinct from radio pulsars 
in their properties, and are believed to be either "magnetars" 
(neutron star s with magnetic fields B > 10 14 G; T hompson & 
Duncan 1996) or exotic accreti ng systems (e.g. v an Paradijs, 
Taam, & van den Heuvel 1995); however, an interpretation 
for the RQNS is less clear. Brazier & Johnston (1999) ar- 
gue that RQNS are energetic young radio pulsars like the Crab 
pulsar, but whose beams do not cross our line of sight. How- 
ever, Vasisht et al. (1997) and Frail (1998) propose that RQNS 
are neutron stars with large initial periods (Pq > 0.5 s) and/or 
high magnetic fields (B ;> 10 14 G) and are thus possibly re- 
lated to the SGRs and AXPs, while Geppert, Page & Zan- 
nias (1999) suggest that they may rather be fast-spinning but 
weakly-magnetized sources. 

One way to distinguish between all these possibilities would 
be to detect pulsations from a RQNS. The period and period 
derivative of th e source could then be use d to infer a surface 
magnetic field ( [Manchester & Taylor 1977[ ), while if the RQNS 
can also be associated with a SNR, an independent age deter- 
mination for the latter can be used to estimate an initial period 



for the neutron star (e.g. Reynolds 1985) 



The RQNS RX J0822-4300 ( petre, Becker, & Winkler 1996| ) 
is near the center of and is almost certain ly associated with 
the young (<5000 yr; |Winkler et al. 1988|; Arendt, D wek, & 
Petre 1991) and nearby (2.2 kpc; feeynoso et al. 1995 ) super- 
nova remnant Puppis A (G260.4-3.3). Recently, Pavlov, Za- 
vlin & Trumper (1999; hereafter PZT99) and Zavlin, Trumper 
& Pavlov (1999; hereafter ZTP99) have analyzed two archival 
ROSAT datasets on RX J0822^1300, separated by 4.6 yr. In 
each dataset they find evidence for weak pulsations, the periods 
of which are slightly different as would be expected for pul- 
sar spin-down. The resulting period, P = 75.5 ms, and period 
derivative, P = 1.49 x 10~ 13 s s , when combined with the age 
of the SNR, imply a dipole magnetic field B = 3.4 x 10 12 G, a 
spin-down luminosity E = 1.4 X 10 37 erg s -1 and an initial pe- 
riod Pq w 55 ms, all of which (despite the radio-quiet nature 
of the source) are properties typical of a young energetic radio 
pulsar associated with a SNR. 

However energetic young pulsars in SNRs have some un- 
mistakable signatures. Every young (<;20 kyr) pulsar located 
within the confines of a SNR powers an observable pulsar wind 
nebula (PWN) — a filled-center synchrotron source resulting 
from the confinement of the relativistic pulsar wind by external 
pressure. Thus a simple test to determine if RX J0822^1300 
is indeed an energetic young pulsar, as argued by Brazier & 
Johnston (1999) and as implied by the detection of pulsations 
by PZT99 and ZTP99, is to see if it h as an associated PW N. At 
radio wavelen gths, existing data (e.g. lArendt et al. 199C|; D ub- 
ner et al. 1991) let us put no useful constraints on the presence 
or absence of a PWN associated with RX J0822-4300. This 
is because these observations were carried out at relatively low 
frequencies (where Puppis A is brightest) and low spatial res- 
olution, resulting in a great deal of confusing emission at the 
position of RX J0822^1300 from both the SNR shell and from 
diffuse internal emission. We have therefore carried out new 
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observations towards RX J0822^300, at higher frequency and 
spatial resolution than previous measurements, aimed at search- 
ing for a radio PWN associated with RX J0822^1300 and thus 
determining whether its properties are consistent with it being a 
young pulsar. Our observations are described in §^|, while in §|| 
we demonstrate the absence of any radio PWN at the position 
of RX J0822-4300 and quantify the consequent upper limits. 
In §|] we argue that this non-detection implies that RX J0822- 
4300 must have properties very different from the young ener- 
getic pulsars which do power observable PWN. 

2. OBSERVATIONS 

Radio observations towards RX J0822^4300 were made with 
the Australia Teles cope Compact Array (ATCA; F rater, Brooks, 
& Whiteoak 1992) in its 0.750D configuration on 1999 July 
24/25. In this configuration, the array contains ten baselines 
in the range 31 m to 719 m, and another five baselines in the 
range 3750 m to 4469 m. Since these two sets of baselines 
cannot be easily combined in a single image, this effectively re- 
sults in two sets of data: one appropriate for imaging extended 
structure on a wide range of scales (a "large scale" image), and 
another sensitive only to a narrow range of spatial scales, but at 
much higher spatial resolution (a "small scale" image). 

Two separate observations were made, each of duration 12 h. 
In the first, data were collected at frequencies of 1.4 and 
2.5 GHz, while in the second, a single observation was made 
centered at 4.8 GHz. At 1.4 and 2.5 GHz the bandwidth was 
128 MHz, while at 4.8 GHz the bandwidth was 256 MHz, in all 
cases divided into 4-MHz channels. 

Observations at 1.4/2.5 GHz consisted of a two-point mosaic 
with mean position centered on RX J0822^300; at 4.8 GHz, 
observations consisted of a single pointing, offset 2.5 arcmin to 
the west of RX J0822^4300 so as to avoid sidelobe contami- 
nation from the nearby bright source PMN J0820-4259. Am- 
plitudes were calibrated by observations of PKS B 1934-638, 
assuming flux densities of 14.9, 11.1 and 5.7 Jy at 1.4, 2.5 and 
4.8 GHz respectively. Instrumental gains and polarization were 
determined using regular observations of PKS B0823-500. 

3. ANALYSIS AND RESULTS 

After standard editing and calibration using the MIRIAD 
package, total intensity images were formed at each frequency, 
using a multi-frequency synthesis approach to both improve the 
// - v coverage and minimize the effects of bandwidth smear- 
ing (|Sault & Wieringa 1994]). At 1.4 and 2.5 GHz, mosaic im- 



ages were formed using maximum entropy deconvolution, both 
pointings being deco nvolved simultaneously (S ault, Staveley- 
Smith, & Brouw 1996). The 4.8-GHz image was deconvolved 
using the CLEAN algorithm. At each frequency, both "large 
scale" and "small scale" images were formed. At 2.5 GHz, the 
two shortest baselines sample emission from the SNR which 
fills the entire field of view of the "large scale" image and pre- 
vents it from being successfully deconvolved; therefore at this 
frequency these two baselines were not used when forming the 
"large scale" image. 

In all cases, deconvolution was constrained to only act on 
specific regions of the image, namely the shell component of 
the SNR (defined usi ng the 0.8 GHz MOST image of B ock, 
Turtle, & Green 1998) and background point sources outside 
the SNR. The resulting model was subtracted from the u- v data 
to produce a dataset which contained visibilities corresponding 

5 The increase in antenna temperature resulting from the flux of the disk is ne 



only to emission from the interior of the SNR. This dataset was 
then imaged, deconvolved, smoothed with a gaussian restoring 
beam, and corrected for the mean primary beam response of the 
ATCA antennas. 

In all six images (1.4, 2.5 and 4.8 GHz; "small scale" and 
"large scale") no radio emission could be seen at or around the 
position of RX J0822-4300; one such image is shown in Fig |l[ 
To quantify these non-detections we performed a series of simu- 
lations, in each of which we modeled the appearance of a PWN 
by using a circular disk of a given surface brightness and ra- 
dius, centered on RX J0822^-300. This simple morphology is a 
reasonable approximation to observed PWN, most of which are 
centered on their associated pulsar with approximately constant 
surface brightness across their extent. In each simulation, the 
Fourier transform of this disk was added to the u-v data from 
which the shell emission had been subtracted, and the imaging 
and deconvolution process was then repeated. 5 For a given di- 
ameter, we increased the brightness of the simulated disk until 
it could clearly be distinguished from the underlying noise (this 
criterion corresponds to a ~ 5a detection for small diameters, 
but is closer to 3a for larger sources). We were thus able to 
quantify the sensitivity of the data to a PWN as a function of 
its size, incorporating effects due to non-Gaussian noise in the 
image, unrelated background sources and the limited range of 
spatial scales sampled by the interferometer. 

The results of these simulations are shown in Fig ||. At each 
frequency, the sensitivity curve consists of four regimes. At 
the smallest scales, each curve is essentially flat, corresponding 
to the sensitivity of the "small scale" image to an unresolved 
source. The curve then increases approximately as S m i„ oc D 2 , 
as ex pected for an extended source (cf. Fig 2 of G aensler et al. 
2000); slight deviations from this relation are due to the effects 
of noise-fluctuations in the data. At a certain scale the sensi- 
tivity of the "large scale" image to an unresolved source be- 
comes better than that of the "small scale" image to a resolved 
source, and the curve becomes flat once more. Finally, at scales 
which are resolved by the "large scale" image, the curve once 
again increases proportional to D 2 . The curve at each frequency 
terminates at the largest scale detectable by the interferometer; 
note that the 2.4 GHz curve ends prematurely because the two 
shortest baselines were not used, as discussed above. 

4. DISCUSSION AND CONCLUSIONS 

To determine whether the limits derived in Fig || are con- 
straining, we need to determine an expected size and flux den- 
sity for a PWN associated with RX J0822^300. A PWN will 
expand until the pulsar wind comes into pressure equilibrium 
with the external pressure. As discussed by Gaensler et al. 
(2000), there are two possible sources for this pressure: either 
the external gas pressure, p gas = nkT, or the ram pressure pro- 
duced by the pulsar's motion, p mm = pV 2 , where n and p are the 
number and mass density respectively of the ambient medium, 
V is the velocity of the pulsar, and T is the temperature of am- 
bient gas. 

By modeling the infrared emission from Puppis A, Arendt, 
Dwek & Petre (1991) derive parameters for the gas interior to 
the SNR (into which the PWN is expanding) of n = 1.3 cm" 3 
and T = 3 - 6 x 1 6 K; similar values are derived from X-ray 
spectroscopy (e.g. Winkler et al. 1981; Berthiaume et al. 1994). 
This corresponds to a pressure ;? gas = 0.05-0.11 nPa which, if 
equated with the pressure E j\-Kr 2 c from the pulsar wind (where 



gligible in all cases. 
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r is the radius of the PWN and E = 1.4 x 10 37 erg s" 1 ), results 
in a PWN of diameter 1 1-16 arcsec at a distance 2.2 kpc. 

For the same number density as used above, and assuming 
the ambient gas to be composed only of atomic hydrogen, we 
find that p mm = 2.2V 3 2 nPa, where V = 10 3 V3 km s" 1 . Balancing 
this pressure with that from the pulsar wind (e.g. G aensler et al. 
2000) results in a PWN of angular diameter ~ W^ 1 arcsec. 

Since the mean velocity of the pulsar population is V3 « 0.38 
(Cordes & Chernoff 1998 ), and in fact in this particular instance 
the offset of RX J0822- 4300 from the dynamical center o f the 
SNR argues that V 3 > 1 ( |Petre, Becker, & Winkler 1996[ ), it is 
likely that p mm > p gas , and that the smaller of the two sizes we 
have just estimated, corresponding to a bow-shock morphology, 
should apply. We note that in such a case, it is still reasonable to 
model the PWN as a circular disk, since in observed bow-shock 
nebulae most of the emission is concentrated close to the head 
of the nebula. In any case, regardless of the dominant source of 
confining pressure, the expected extent of a PWN powered by 
RX J0822^4300 is small. Although we believe the sizes derived 
above to be robust, we conservatively adopt a maximum angu- 
lar size for any PWN of 30 arcsec to take into account possible 
uncertainties in V, n, T, or the distance to the source. From 
Fig ||, it can be seen that at all three frequencies, the upper limit 
on the flux density for such a source is ^7 mJy. 

Assuming a typical PWN spectral index of a = -0.3 (JS V oc 
v a ), an upper limit of 7 mJy at 1.4 GHz corresponds to a 
broad-band radio luminosity (integrated between 10 MHz and 
100 GHz) of L R = 2 x 10 30 erg s _1 . Defining e = L R /E to be 
the ratio between a PWN's broad-band radio luminosity and its 
spin-down luminosity, we find that for E = 1.4 x 10 37 erg s -1 
as reported by PZT99 and ZTP99, we can derive an upper limit 
e < 10~ 7 . This is a more stringent limit on e than has been de- 
rived for almost any other pulsar (cf. Frail & Scharringhausen 
1997; paensler et al. 2000| )! In particular, this upper limit 
is sharply at odds with the values of e seen for other young 
(<; 20 kyr) pulsars, almost all of which produce radio PWN or 
have upper l imits consistent with e > 10" 4 (F rail & Scharring- 
hausen 1997; paensler et al. 2000| ). The glaring exception to 
this is PSR B1509-58 in the SNR G320.4-1.2 (MSH 15-52), 
which powers an X-ray PWN but for w hich no radio PWN has 
yet been detected ( paensler et al. 1999 ). However, this can be 



understood in terms of the low ambient density (« < 0.01 cm" J ), 
which results in severe ad iabatic losses and a consequently un- 
derluminous radio PWN ( Bhattacharya 1990 ). This condition 
is not satisfied for RX J0822^-300, and so cannot be consid- 
ered as a possible explanation for the non-detection of a PWN. 6 
We thus find that any PWN in Puppis A has a radio luminos- 



ity three orders of magnitude fainter than expected for the spin 
parameters derived by PZT99 and ZTP99. 

Nevertheless, if we assume, as Brazier & Johnston (1999) 
have argued, that RX J0822-4300 is a rotation-powered pul- 
sar, what spin parameters can we infer for it? If we require that 
e > 10~ 4 as seen for other young pulsars, the maximum value of 
E = 4n 2 IP/P 3 which is consistent with our non-detection of a 
radio PWN is ~ 10 33 erg s" 1 . Meanwhile, it is unlikely that the 
characteristic age, r = P/2P, of the pulsar is more than 50 kyr, 
~10 times the true age of the system. These upper limits on E 
and t correspond to lower limits P > 3.5 s, P > 1.1 X 10~ 12 s s" 1 
and B>6.4xl0 13 G, parameters which are very similar to 
those seen for the SGRs/AXPs (K aspi, Chakrabarty, & Stein- 
berger 1999) and for the young radi o pulsar PSR J1814- 1 744 
dPivovaroff, Kaspi, & Camilo 2000t |Camilo et al. 2000| ), but 
quite different than those of other young pulsars in SNRs, for 
which typically P < 0.2 s, E > 10 36 erg s" 1 and B w 10 12 G. 

Whether RX J0822^4300 indeed has a long initial period and 
high magnetic field, or has some other properties such that it 
does not produce a detectable radio nebula or radio pulsations, 
the lack of a PWN around this source (and around other RQNS 
such as IE 1207.4-5209 in the SNR G296.5+10.0; M ereghetti, 
Bignami, & Caraveo 1996; piacani et al. 2000| )~argues that 
at least some RQNS have drastically different properties from 
young radio pulsars. 

Brazier & Johnston (1999) list six RQNS which are younger 
than 20 kyr and nearer than 3.5 kpc. Excluding two RQNS 
from their list which do power PWN and thus may well be ra- 
dio pulsars beaming away from us, but including the recently- 
discovered RQNS in the young and nearby SNR Cassiopeia A 
( [Tananbaum 1999j ; |Pavlov et al. 2000| ; phakrabarty et al. 2000| ), 
this implies a Galactic birth-rate for such sources of at least 
once every ^200 years, compar able to or even in excess of the 
birth-rate for radio pulsars (e.g. Lyne et al. 1998 ). Radio-quiet 
neutron stars thus point to the possibility that pulsars like the 
Crab are not the most common manifestation of neutron star. 
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FIG. 1. — "Large scale" ATCA image of the region surrounding RX J0822-4300, at a frequency of 2.5 GHz Emission frnm the shpll nf Puppjs A and from 
bright point sources exterior to the SNR has been subtracted. The circle marks the position of RX J0822^4300 ( Petre, Becker, & Winkler 199C ) — the positional 
uncertainty is half the radius of the circle. The greyscale ranges from 0.0 to 2.4 mjy beam -1 at a resolution of 27 X 17" (indicated at lower right); the rms noise is 
0.2 mjy beam -1 . 
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FIG. 2. — Minimum detectable flux density of a pulsar wind nebula at the position of RX J()822^13()0 at 1 .4, 2.5 and 4.8 GHz, as a function of size. 



